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^ ! Abstract: A bosonic reahzation of the SU{2) algebra and of its vector representation is 

^1 constructed, and an effective low-energy description of the Bose-Hubbard model in the form 

c/3 ; of anisotropic theory of quantum rotors is proposed and discussed. A possibility of a normal 

-— . ■ zero-temperature bosonic phase with neither crystalline nor superfiuid order around the tip 

of the checkerboard-solid lobe at half-integer fillings is examined. 
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The ground state of interacting bosons is typically a superfluid, unless a sufficiently long- 
ranged interparticle interaction or a commensurate periodic potential is present, in which 
case bosons can form a crystalline solid even at T = 0. The simplest description of the 
interacting bosonic systems is provided by the Bose-Hubbard model and its variants, which 
depending on range and strengths of interactions and on density of particles exhibit uniform 
Mott insulating, charge-density- wave (solid), superfluid, and mixed supersolid ground states 
im. Influence of quantum fluctuations on the phase diagram of quantum bosons has been 
somewhat less well understood, and recently a possibility of a normal phase of bosons at 
T = has been raised in the context of experiments on superconductor-insulator transition 
in thin films 0. The normal phase would be without either solid or superfluid order, and, 
presumably, would show the conducting properties of a metal. Leggett has supplied argu- 
ments against such a ground state for bosons in the continuum some time ago 0, but on a 
lattice the situation is less clear. The issue is not only of theoretical interest, as there exist 
many experimental realizations of lattice bosonic systems, including granular superconduc- 
tors 0, Josephson junction arrays 0, or, since recently, cold atoms in optical lattices [||. 
Recent numerical calculations which focused on this question 0], ||^ in D = 1 gave results 
that seem to depend on the details of the model, and no coherent picture has yet emerged. 

It is well known, on the other hand, that for infinite on-site repulsion, which limits the 
boson occupation numbers at each site to zero and one, bosonic system becomes equivalent 
to the system of spins-1/2. Das and Doniach p| argued that in the opposite limit of a 
large number of particles per site, the equivalent quantum-phase (or Josephson junction 
array) model can be written in terms of the generators of the Euclidean group E2, which 
again is intimately related to, in fact just a simple contraction of, the group 5*0(3). It 
appears natural to ask if there exists a more general connection between the rotational 
symmetry and the lattice bosons. One such connection is established in this paper, where 
I display a bosonic realization of the SU{2) Lie algebra and show that the two-component 
superfluid order parameter and the order parameter for the checkerboard solid phase can be 
combined into this algebra's three dimensional vector representation. Similar operators were 
constructed previously for the electronic Hubbard model 0, and applied to its negative-U 
version. Here, it is proposed that near half-integer fillings the Bose-Hubbard model with on- 
site and nearest-neighbor repulsion may posses an approximate SU{2) symmetry, and that 
the effective low-energy theory near these points has the form of anisotropic Hamiltonian for 
quantum rotors. At half-integer ffilings, where the mean-field theory predicts a ffist-order 
(spin-flop) transition between solid and superfluid phases, the effect of quantum fluctuations 
is the strongest, and right at the boundary it can produce a disordered phase, with only 
short-range crystalline or superfluid order. Since such a phase has a finite gap and therefore 
is stable to small perturbations, it will occupy a finite region of the phase diagram. In 
this case the direct solid-to-superfluid transition is replaced by two continuous (of the Ising 
and XY variety) transitions between solid and normal, and normal and superfluid phases. 
In the example of hard-core bosons with nearest-neighbor and weak next-nearest-neighbor 
interactions, the effective rotor theory suggests that the intermediate normal phase is most 
likely absent in dimensions D > 1. 1 speculate that in D = 1 the normal phase may be 



observable at larger filling factors or for weaker on-site repulsion. 

Consider a system of bosons on a quadratic lattice in Z)- dimensions and at T = 0. It is 
straightforward to check that the operators: 

p 

L- = '^T.bpb-,^Q^ (2) 

p 

L3 = l{T.bpP + ^l (3) 

p 

form the SU{2) Lie algebra: [L3,L±] = ±L±, [L+,L_] = 2L3. Here [bi,bj] = 5ij are the 

standard Bose operators, 6, = J^p^wiW' Ri)bp, Q = {Tf/a, n/a....'7i/a) and the sums over mo- 
menta run over the first Brillouin zone. Nl is the total number of lattice sites. Components 
of the three dimensional vector representation of SU{2) should satisfy the commutation re- 
lations: [La, Up] = ieaiB'yn^, where a = 1,2,3, and L± = Li ±iL2. A possible SU{2) vector 
multiplet is then: 

^i = ^E(^M + W, (4) 

i 

n2 = ^Y.(b\b\-hh), (5) 



«3 



^-Y^dQH%. (6) 



The third component is proportional to the amplitude of the charge- density- wave with the 
period of two lattice spacing in every direction, and represents the order parameter for the 
checkerboard-solid phase. It is interesting that in order to satisfy the SU{2) commutation 
relations it is necessary for rii and ^2 to be bilinear in Bose operators. The finite expectation 
value of ni or 77-2 thus signals the long-range order in the composite operator bb, not b itself. If 
there is no attractive interaction in the microscopic Hamiltonian that could lead to pairing 
of bosons it seems plausible to assume that the composite operator bb can condense only 
when (6) 7^ 0. I will take this to be the case and hereafter make no difference between the 
two forms of superfiuid order. 

In general, the Bose-Hubbard Hamiltonian 

H = -tY, blbj + 2 E ^ij^if^j - l^Yl^i, (7) 

of course commutes with L3, but not with Li and L2, and does not posses the exact SU{2) 
symmetry. However, the numerical calculations and the mean-field theory |]l| show that right 
at half-integer fillings the typical system described by the Bose-Hubbard model exhibits solid 



and superfluid phases at T = 0, which brake Ising and U{1) symmetry, respectively. The 
discrete Z2 and the f/(l) symmetries can be embedded into a larger SU{2) group, and the 
existence of the bosonic version of the SU(2) algebra and of its vector multiplet allows one 
to construct a phenomenological theory for the Bose-Hubbard model that would apply at 
low energies and provide a unified description of solid and superfluid order. In continuum 
notation, the effective Hamiltonian density assumes the form: 

Heff = -E,i[n] + -i—Ll + PaiViU^f) - KL^, (8) 

where n and L are dependent on the D-dimensional coordinate x and the imaginary time r, 
and the Hamiltonian is obtained by integration over space and imaginary time, as usual. L 
can be thought of as the angular momentum. We allow that the SU{2) symmetry is broken to 
f/(l) X Z2 both implicitly, by xi,2 = X± 7^ Xa ^"^^ Pi, 2 = P±i^ Ps, and explicitly, by including 
the source fields h and the classical energy Ed, which favor a particular orientation of the 
vector n. I assume the simplest form Ed = gn\. Also, the fluctuations of the order parameter 
amplitude will be taken to represent high-energy modes, and only the orientation of n then 
remains at low energies. Thus n^ = 1, and L ■ n = 0. 

Since the rotor Hamiltonian (8) is quadratic in L one may integrate the generators out 
to obtain the effective theory for the order parameters only. This requires some care because 
of the orthogonality between L and n. I introduce a "gauge-field" as L = n x A, and add 
a " gauge- fixing" term to the Hgff as {l/2uj){'n ■ A)"^. The integration over A can then be 
performed, and it results in the Lagrangian density 

Leff = -Ed[n] + -{paiViU^f + x±(^ -inx hf + (9) 

X± 5-(n X n — in X n X h)'i), 

1 — vnj_ 

with f = 1 — (xx/xs); independent of the gauge-fixing parameter u. Similar form has been 
studied by Zhang in context of his SO (5) theory of high- Tc superconductivity flOl, but even 



in the classical approximation n = Vn = the Eq. 9 differs from his expression ||rT|. Only 
if the symmetry breaking is weak, i. e. f << 1, the classical ground states of the system are 
determined by minimization of a simpler expression equivalent to Zhang's: 

1 h"^ ^ 

Ld ~ -gnl - —n\{xLnl + Xs^l), (10) 

where xi. = X±, Xz = X±(l + '^); and /i = (0, 0, /i). At /i = there is a first-order (spin-flop) 
transition from solid to superfluid as g changes sign from negative to positive. At finite h 
and if Xs > X-L? the same remains true only at negative g, while if xs < X± there exists an 
intermediate supersolid phase in between. 

The effect of quantum fluctuations on the classical solutions is expected to be the 
strongest right at the point g = h = 0, where the classical energies of solid and super- 
fluid phases become degenerate. After rescaling the imaginary time as {p±/x±y^'^T ~^ t the 



fluctuation part of tfie Lagrangian becomes: 

Leff = 7^7^^((Vin,)2 + {ky + uiViTisy + vin X fi)l), (11) 

wfiere u = {pz/ Pl) — 1- It is tlie non-linear sigma model in D + 1-dimensions, at an effective 
temperature Tejj = l/(p±X±)^''^) with two new terms representing the remaining implicit 
symmetry breaking. To determine the phase of the system at the point /i = I perform 
the renormalization group transformation on the model, near D = 1 and perturbatively in 
Tg//, u and v |jl2|]. The integration over the degrees of freedom within the momentum shell 



(A/s, A) changes the couplings of the model according to the differential equations: 

dW) ^ ^ ^ '" 241 + M) ^^ff ^^^> 

^^-^ = --{3.,(m-2MT.,„ (13) 

'^''' -2M-4^, (14) 



rfln(.) '^' Al + \9\y 
where the general number of components m is retained, m = 3 being the case in point. 6{g) 
is the step function, and I redefined gT^ff/A'^ -^ g. Weak implicit symmetry breaking is, 
expectedly |TB[, |T^, irrelevant, and in D > 1 aX g = the system has a phase transition 



at the critical effective temperature T*^^ = 2n{D - 1) + 0{{D - 1)^). For T^jf < T*^f 
the fluctuations are too weak to change the mean- field result, and there is indeed a direct 
first-order transition between solid and superfluid at (7 = 0. In the opposite case quantum 
fluctuations qualitatively change the phase diagram. Since the disordered phase, which then 
exists right at (/ = /i = 0, has a finite correlation length and thus a finite gap, it can not 
collapse at infinitesimal g or h. Consequently, the disordered phase spreads from the point 
(yf = /i = to a finite region around it, and gives way to the solid or the superfluid only at 
some finite values of g or h. This also becomes evident from the above recursion relations 
at finite g, which reduce to those in ref. 13 when (7 > 0. The resulting flow diagram for 
h = u = V = Q is depicted for m = 3 and Z) > 1 at Figure 1. By increasing g at fixed 
Teff > T*jrj one intersects both solid-to-normal and normal-to-superfluid phase boundaries, 
and the normal phase persists for a finite interval of g. 

Whether or not will a particular example of the Bose-Hubbard model exhibit the 
intermediate normal phase depends on the effective temperature at the symmetric point 
g = h = 0. In general, there exists no simple way to calculate the effective coupling con- 
stants in terms of the microscopic interactions. For the hard-core bosons with infinite on-site 
repulsion, however, the connection between microscopic and effective couplings can be made 
more explicit, due to the well known analogy with the spin-1/2 system. Taking Vij = V for 
{i,j) nearest neighbors, and Vij = V << V for the next-nearest neighbors, one can derive 
the effective rotor Hamiltonian in Eq. 8 in large spin-S approximation |TB], and eventually 
take S = 1/2. This way one finds: 

Ea = -S\2ztnl + {zV - z'V')nl), (15) 



h = fi-^{zV-z'V'), (16) 

X-,} = 2zt-S-'E,i, (17) 

X^' = zV + z'V'-S~^E,i, (18) 

P± = S\ (19) 

P3 = y(l^-n- (20) 

z and z' are the number of nearest and next-nearest neighbors on the lattice, respectively. 
For V = 0, inserting these expressions into the Eq. 9, the classical approximation yields the 
first-order transition from solid to superfluid at: 



tc = -^Kv-f^) (21) 

for 5* = 1/2, in exact agreement with the result of the mean-field approximation applied 
directly to the Bose-Hubbard Hamiltonian, in two dimensions |]T[]. In this case at half- 
filling (at /i = 2V) and at t = V/2 the system has the full rotational symmetry, and in its 
immediate vicinity Xs = X±j so that there is no supersolid phase. When small y > is 
turned on, exactly at half-filling (now at /x = {zV — z'V')/2) the solid-to-superfluid transition 
at 2t = V — z'V /z remains first-order, but since xs = X-l(1~^'^V^^) < X±j away from half- 
filling there immediately appears a thin sliver of the supersolid phase in between. Again, 
this is precisely what the direct mean-field treatment of the Bose-Hubbard model in the 
hard-core limit obtains, and our effective theory provides a very simple explanation for the 
known mean-field phase diagram. 

Given the good agreement at the mean-field level we may attempt to estimate the 
effective temperature for hard-core bosons at the transition point at half-filling. Near D = 1 
where the effective temperature is dimensionless, for S = 1/2 I find Te// = 2\/2z, which 
suggests that in D > 1 hard-core bosons with weak next-nearest-neighbor repulsion are 
probably too "cold" to have the intermediate disordered phase. For example, in D = 2 
Teff ~ 4^/2 while T*fr ^ 27r, according to this very crude estimate. The case of D = 1 is 
special, since the presence of additional topological terms in the effective action makes the 



disordered phase at the solid-superfluid boundary gapless for half- integer spins |T5[ . Instead, 
in D = 1 one finds the continuous Kosterlitz-Thouless solid-superfluid transition right at 
half-filling |1[], [Q. If the constraint on double occupancy is relaxed the Bose system is no 
longer exactly equivalent to the spin-1/2 Heisenberg Hamiltonian, but one would expect that 
the form of the low-energy theory should remain the same. For smaller on-site repulsion and 
larger half- integer fillings the multiple occupancy of the same site should become frequent, 
and it is conceivable that the rotor Hamiltonian then may become appropriate even in 
D = 1. This would suggest that the under these conditions the intermediate normal phase 
may appear at the tip of the solid-superfluid lobe, in qualitative agreement with the Monte 
Carlo calculations on the quantum-phase model of ref. 7 ||16| . 



The normal phase of bosons discussed here differs in detail from the one proposed by 
Das and Doniach 0, which appears at integer filling in the quantum phase model for large 
enough nearest-neighbor repulsion. Both phases are the results of quantum disordering of 
the bosonic ground state, however, and should show similar physical properties. Studies of 
the conducting properties of this phase and of the effects of disorder [^ are left for future 
research. 

In conclusion, I presented a bosonic representation of the SU{2) Lie algebra, and showed 
that the superfluid and solid order parameters together form its vector representation. It was 
argued that the low-energy effective theory of the Bose-Hubbard model near half-integer fill- 
ings has the form of anisotropic Hamiltonian for quantum rotors. In classical approximation, 
the phase diagram derived from the effective theory agrees well with the previous studies. 
A possibility of a quantum disordered, or a normal, phase between solid and superfluid is 
pointed out, and the conditions for its appearance are briefly discussed. 

This work has been supported by NSERC of Canada. 



FIGURE CAPTIONS: 



Figure 1: Schematic flow and phase diagram of the effective rotor theory at half-integer 
filling h = 0. At high effective temperature T^jf, with tuning of the parameter g the system 
undergoes the solid-to-normal phase transition governed by the Ising critical point in D + 1, 
followed by the normal-to-superfiuid transition in the XY universality class. At low Tg/j the 
trivial fixed point governs the direct first-order solid-to-superfiuid transition. 
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